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This work describes a new microfluidic device developed for rapid screening of solubility diagrams. 
In several parallel channels, hundreds of nanoliter- volume droplets of a given solution are first stored 
with a gradual variation in the solute concentration. Then, the application of a temperature gradient 
along these channels enables us to read directly and quantitatively phase diagrams, concentration 
vs. temperature. We show, using a solution of adipic acid, that we can measure ten points of the 
solubility curve in less than 1 hr and with only 250 /iL of solution. 



I. INTRODUCTION 

Chemistry, biology, and pharmacology, are facing al- 
ways more complex systems depending on multiple pa- 
rameters. Therefore their complete investigations take 
time and require significant amounts of products. In this 
context, robotic fluidic workstations have already met a 
great success and proved their efficiency for instance in 
the genome sequencing and analysis However, these 
instruments remain very expensive, need important la- 
bor, and the volumes involved (< mL) are still too large 
for some specific applications (e.g. proteomics) 0,0] ■ 

Nowadays, other high throughput techniques based on 
microfluidics 0, Q can offer suitable alternative solutions 
for the development of rapid screening tools. Microfluidic 
devices are now largely used in biological and chemical 
fields for multiple applications like molecular sepa- 
rations and cells sorting Q, polymerase chain reaction 
[1, Hi, rapid micromixing and analysis of chemical reac- 
tions [l(J, [ll|, EH, Moreover, the development of 
microvalves and micromixers has made possible the pro- 
duction of highly integrated systems which can be used 
to address individually hundreds of reaction chambers 
These devices are well adapted to carry out high 
throughput screening of phase diagrams, particularly in 
the case of protein crystallization investigation. However, 
their fabrication and multiplexing are still complicated. 
Another possible strategy is the use of droplets playing 
the role of nanoliter-sized reaction compartments. These 
droplets can be produced in specific microfluidic geome- 
tries and their volume and chemical composition can 
be fixed in a controlled way. In addition, they also allow 
a rapid mixing of the different compounds, prevent from 
any hydrodynamic dispersion and cross contamination, 
and can be stored in microchannels (see Ref. [IB] and 
references therein). Such a strategy has already proved 
to be useful for crystallization studies: eji. screening 
of protein crystallization conditions pj], El, or crystal 
nucleation kinetics measurements 

Figure Q] summarizes the main insights of our work. 
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We have engineered a new microfluidic chip that allow 
a direct and quantitative reading of two-dimensional di- 
agrams. Hundreds of nanoliter-sized droplets of differ- 
ent chemical compositions can be stored in parallel mi- 
crochannels, and a temperature gradient applied along 
these channels enables us to obtain a two-dimensional 
array of droplets of different concentrations and tempera- 
tures. For solubility diagram screening, droplets contain- 
ing a given solute are first stored with a gradual variation 
of concentration. Then, crystallization in the droplets is 
induced by cooling, and finally, the application of an ad- 
equat temperature gradient dissolves crystals in droplets 
whose temperature is higher than their solubility tem- 
perature. As a result, we directly read the limit between 
droplets with and without crystals as shown on Fig.QJc), 
which gives the solubility temperatures of the solution at 
the different concentrations. 

In the materials and methods section, we describe the 
microfluidic device, the method used to store the droplets 
in the channels, and the temperature control setup. We 
also characterize the concentration and temperature gra- 
dients. In the last section, we present an experimental 
protocol to measure quantitatively solubility diagrams 
using this device. We demonstrate its efficiency by mea- 
suring with only 250 /iL of solution, the solubility curve 
of an organic compound. 



II. MATERIALS AND METHODS 

A. Microfabrication 

The microfluidic device is fabricated in 
poly(dimethylsiloxane) (PDMS) by using soft- 
lithographic techniques [2fl]. PDMS (Silicone Elastomer 
Base, Sylgard 184; Dow Corning) is molded on master 
fabricated on a silicon wafer (3-Inch-Si- Wafer; Siegert 
Consulting e.k.) using a negative photoresist (SU-8 2100; 
MicroChem). To make molds of 500 /im height, we spin 
successively two 250 /im thick SU-8 layers on the wafer. 
After each spincoating process, the wafer is soft-baked 
(10 min/65°C and 60 min/95°C). Photolithography is 
used to define negative images of the microchannels. 
Eventually, the wafer is hard-baked (25 min/95°C) and 
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FIG. 1: (a) Design of the microfluidic device (channels width 
500 ^m). Silicone oil is injected in inlet 1 and aqueous so- 
lutions in inlets 2 and 3. The two dotted areas indicate the 
positions of the two Peltier modules used to apply tempera- 
ture gradients VT. The three lines of dots mark the positions 
of temperature measurements, (b) Picture of the microfluidic 
chip made of PDMS sealed with a glass slide (76x52 mm 2 ) to 
improve clarity. Droplets containing a colored dye at different 
concentrations are stored in the ten parallel channels, (c) Ex- 
ample of direct reading of a solubility diagram. The droplets 
contain an organic solute. The dotted line bounding droplets 
containing crystals give an estimation of the solubility limit 
(see section Results for details). 



developed (SU-8 Developer; MicroChem). A mixture 
10:1 of PDMS is molded on the SU-8 master described 
above (65°C/60 min). The crossed linked PDMS layer 



is then peeled off the mold and holes for the inlets 
and outlets (1/32 and 1/16 in. o.d.) are punched into 
the material. Then, the PDMS surface and a clean 
silicon wafer surface (3-Inch-Si-Wafer; 500 /tm; Siegert 
Consulting e.k.) are activated for 2 min in a UV ozone 
apparatus (UVO Cleaner, Model 144AX; Jelight) and 
brought together. Finally, the device is placed at 65° C 
for 2 hr to improve the sealing. 



B. Droplet storage protocol 

The device, presented on Fig. [Ha), is composed of 
three inlets and ten outlets located at the extremities 
of channels Ci to Cio- As shown on Fig. [ljb), each out- 
let is connected to a ~ 20 cm long rigid tubing (FEP 
1/16 in.) ended with a piece of soft PVC tubing (Nal- 
gene, « 5 cm long) inserted in an automated pinch elec- 
trovalve (105S-01059P; Asco Joucomatic). Thanks to 
this system, each outlet can be independently closed or 
opened by pinching or not the corresponding PVC tub- 
ing. However, the pinching out of a tube leads in a liquid 
displacement. To minimize the subsequent liquid distur- 
bance in the microchannels, the electrovalves are placed 
close to the rigid ones, and the hydrodynamic resistance 
after the electrovalves is kept as weak as possible using 
large tubing. 

Silicone oil (500 cSt; Rhodorsil) is injected in inlet 1 at 
constant flow rate Q\ « 3 mL hr -1 , and aqueous phases 
are injected at flow rates Q2 and Q3 ranging from to 
about 1 mL hr -1 , in inlets 2 and 3 respectively. All 
liquids are injected with syringe pumps (PHD 2000 infu- 
sion; Harvard Apparatus). At the intersection between 
the oil and the aqueous streams, monodisperse droplets 
of the aqueous phase in oil are continuously produced 
[1H. Both the droplet volume (about 100-300 nL) and 
their production frequency (typically between one to ten 
droplets per second) can be tuned by the ratio of oil to 
aqueous phase flow rates. The droplet composition is 
monitored by the ratio Q2/Qi- 

Thanks to the possible opening and closing of each out- 
let, we can store droplets of given aqueous compositions 
in the different storage channels C4. Several steps are nec- 
essary to perform such a filling. First, all the channels 
are initially filled with silicone oil. Secondly, the outlet of 
channel ci is opened and all the others are closed. In this 
configuration, all the droplets of a given composition flow 
through ci- Finally, once the flow is stable, the outlet of 
Ci is suddenly closed and simultaneously, the outlet of 
channel C2 is opened. All the droplets previously present 
in ci stay immobilized whereas the other droplets, whose 
composition can be changed, flow through C2. Succes- 
sively, in the same way, we can store droplets of various 
compositions in all channels c,. 
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C. Chemical composition control 

For solubility investigations, the control of the con- 
centrations in the droplets is crucial. However, because 
of PDMS elasticity and syringe pumps precision, an in- 
accuracy in droplet concentration remains. To estimate 
this error, we have performed investigations with a con- 
focal Raman microscope (HR800 Horiba; Jobin-Yvon). 
A 50 x microscope objective was used for focusing a 
532 nm wavelength laser beam in the droplets, and for 
collecting Raman scattered light, subsequently dispersed 
with a grating of 600 lines per millimeter. To minimize 
the out-of-focus background signals, we fixed the con- 
focal pinhole at 500 fim. Experiments were performed 
on droplets made of two initial aqueous solutions of 
K 4 Fe(CN) 6 (0.5 M) and K 3 Fe(CN) 6 (0.5 M) injected in 
inlets 2 and 3 respectively. These two compounds display 
strong and distinct Raman signals j22j. 

Figure [2] shows three Raman spectra measured 
in droplets containing different concentration ratios 
Rc=[K 4 Fe(CN) 6 ]/[K 3 Fe(CN) 6 ]. The two bands centered 
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FIG. 2: Raman spectra of droplets containing different con- 
centration ratios Rc of potassium ferrocyanide K4Fe(CN)6 
and potassium ferricyanide K3Fe(CN)6- (a) Rc = (b) 
Rc = 1 (c) Rc = 9. 

at 2060 and 2095 cm' 1 correspond to K 3 Fe(CN) 6 and 
the one at 2136 cm -1 corresponds to K 4 Fe(CN) 6 . The 
concentration of each compound can be probed from the 
area under their specific Raman bands by: 

At = KidtV , (1) 

where is the area under the Raman band of the com- 
pound i, Ci its concentration, Ki a specific constant, t 



the acquisition time, and V the analysis volume. As a 
consequence, the ratio Ra of the Raman bands areas of 
K 4 Fe(CN)6 and K 3 Fe(CN)6 is proportional to the con- 
centrations ratio Rc, and does not depend on the acqui- 
sition parameters. 

In order to optimize the filling protocol, we first use Ra- 
man microscopy to follow the kinetics of the concentra- 
tion stabilization in the droplets after a sudden change in 
the aqueous phases flow rates. Indeed, due to the PDMS 
elasticity and the injection system (syringe pumps), the 
finite response time of the device does not allow instan- 
taneous change of the concentrations. To estimate this 
response time, we have performed the following experi- 
ment: for t < s, Q2 — and Q 3 = 500 fj.L hr _1 , and 
for t > s, Q 2 = Q3 = 250 hr" 1 . Droplets first 
flow through channel ci which is closed after 30 s. Then, 
droplets are stored in five other channels after 1, 2, 4, 
6, and 10 min. Thus, Raman spectra obtained from the 
droplets in the different channels enable us to follow the 
evolution of Ra as a function of time after the flow rates 
change. Figure [3](a) shows it reaches almost a constant 
value after 60 s meaning the concentrations become sta- 
ble after this time. Such measurements illustrate that 
20 min long protocols are efficient to store droplets of 
desired compositions in the ten channels (« 2 min per 
channel). 

A second series of experiments was performed to es- 
timate and characterize the concentration gradient we 
can apply in the device. The storage channels are filled 
with droplets of different concentrations in K 3 Fe(CN)g 
and K 4 Fe(CN) 6 set from the flow rates. In each channel 
Ci, to reach a stable droplets composition, we maintain 
the flow for 90 s before closing the outlet to store them 
[see Fig. HJa)]. By measuring the Raman spectra of the 
droplets composition in the different channels, we ob- 
tain the ratio Ra as a function of the theoretical ratio 
Rc = Q^IQz- The error bar corresponds to the stan- 
dard deviation of the measurements performed on the 
droplets in a given channel. As can be seen on Fig. 03(b), 
a linear relationship between Ra and R c is observed as 
expected. Deviations of a few percents around the linear 
law are probably due to the Raman measurements un- 
certainties, to the accuracy of the injection system, and 
also to the PDMS elasticity. 

These Raman measurements demonstrate that with 
the developed protocol, we are able to store hundreds 
of droplets in ten channels in about 20 min, and con- 
suming less than a few hundreds of [iL of solution. We 
believe that more rigid and smaller microdevices com- 
bined with even more reactive injection system would 
decrease signicantly the amount of liquids used when fill- 
ing the channels. Other strategies involving for instance 
droplet generation thanks to integrated microvalves [13], 
may also proved to be useful to decrease the required 
volumes of solution. 
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FIG. 3: (a) Evolution of the ratio Ra in the droplets after a 
sudden change of the aqueous solutions flow rates Q2 and Q3. 
Before t = s, Q2 = and Q 3 = 500 fiL hr _1 . For t > s, 
Q2 ~ Qs = 250 /iL hr _1 . Between 5 and 30 s, Ra are obtained 
from three single droplets in channel ci . After t = 60 s, each 
point is a mean value calculated on several droplets in a given 
channel, (b) Concentration ratio Ra in droplets as a function 
of the concentration ratio Rq determined from the aqueous 
solutions flow rates. The dotted line corresponds to the linear 
fit of the data. 



a data acquisition instrument (USB-9161; National In- 
struments) and Lab View software. Figure [Ufa) shows 
we are able to apply easily temperature gradients up to 
45°C on 5 cm. To estimate the temperature at any po- 
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D. Temperature control 

The temperature field of the chip is controlled with two 
Peltier modules (30x30x3.3 mm 3 ; CP1.4-71-06L; Mel- 
cor) placed underneath the wafer at positions marked by 
the two dotted areas on Fig. QJa). Since the two Peltier 
modules are independant, we can heat or cool the de- 
vice, and also apply important temperature gradients. 
We use a silicon wafer as chip support to optimize ther- 
mal transfers and thus to create regular temperature gra- 
dients along the storage channels. Thin thermocouples 
(type K, 76 /xm o.d., 5SRTC-TTKI-40-1M; Omega) mea- 
sure the temperature of the device along three series of 
positions parallel to the storage channels. The first series 
is placed above ci, the second one between C5 and eg, and 
the third one below cio [see Fig.QJa)]. To reach the maxi- 
mal precision on liquid temperature measurements inside 
the channels, the thermocouples are inserted in holes pre- 
viously punched through the PDMS layer and filled with 
silicone oil. Thermocouples signals are processed with 



FIG. 4: Temperature profiles of the chip for a given tem- 
perature gradient (a) Temperatures measured along the stor- 
age channels with thermocouples inserted through the PDMS 
layer at different positions shown on Fig. [l]Ja). (A) measure- 
ments series above channel ci; (■) series between cs and C6; 
(T) series below cio. (b) Interpolated temperature profile of 
the chip. 



sitions along the storage channels, we perform a longi- 
tudinal and transverse linear interpolation of the three 
series of measurements. The final profile obtained after 
such interpolation is depicted on Fig.[HJb). Note that the 
temperature is not perfectly homogeneous transversely 
to the storage channels. This is due to the size of the 
Peltier module as compared to the size of the droplet 
storage area: smaller storage area, or larger Peltier mod- 
ules, would give homogeneous temperature profiles along 
the transverse direction of the channels. 
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III. RESULTS 



In the previous section we have shown that our mi- 
crodevice allows us to build a two-dimensional array of 
droplets with both concentration and temperature gra- 
dients. We now present an application of this chip by 
measuring the solubility curve of an organic solute. 

Such measurements are carried out with an adipic acid 
solution previously prepared in a beaker. It is made of 
10.14 g of adipic acid (99%; Aldrich) in 50.66 g of deion- 
ized water. The solubility temperature of the solution 
is 63° C. To avoid any crystallization before the droplets 
formation, the syringe containing the solution and the 
corresponding tubing are heated at about 65° C with two 
flexible heaters (Minco) controlled with temperature con- 
trollers (Minco). A stereo microscope (SZX12; Olympus) 
with an objective (DF PLFL 0.5 x PF; Olympus) enables 
us to observe the device during the solubility study. 

We inject the adipic acid solution in inlet 2 and deion- 
ized water in inlet 3. By changing the flow rates ratio 
we fill the storage channels with droplets whose concen- 
tration in adipic acid varies from 20 g / 100 g of water 
in channel Ci down to 6 g / 100 g of water in Cio- The 
massic concentration C in the droplets is calculated ac- 
cording to: 

C° 

° = 1 + (1 + C°)Q 3 /Q 2 (2) 

where C° is the massic concentration of the initial adipic 
acid solution, Q 2 and Qs the respective flow rates of the 
solution and water (we checked that density variations in- 
duced by the presence of adipic acid are negligible) . The 
microfluidic chip is kept at about 65° C using the Peltier 
modules to avoid any crystallization during the droplet 
storage. Before stopping the droplets in a channel, we 
maintain it open for 90 s for flow stabilization. In these 
conditions, the total filling of the ten channels is reached 
in less than 20 min and only 250 /uL of solution are spent. 

After the droplet storage, crystallization is induced by 
cooling. Note that the mean time of crystal nucleation 
is inversely proportional to the reactor volume. Indeed, 
the nucleation frequency is given by 1/JV where J is the 
nucleation rate that does no depend on the volume V 
of the reactor (see Refs. [IsL [HI, HH, HB| and references 
therein). Crystal nucleation in a droplet of 100 nL is 
thus 10 4 times longer than in a vial of 1 mL. To reduce 
such long induction time, we apply a strong cooling to 
increase significantly the supersaturation. In our case, 
down to w — 5°C, crystals appear in all the droplets after 
a few minutes. 

To obtain the solubility curve directly on the chip, we 
then apply a temperature gradient between 32 and 65° C 
after the crystallization step. Crystals dissolve in all the 
droplets whose temperature is higher than their solubil- 
ity temperature. In the other droplets, crystals are partly 
solubilized but still exist (the equilibrium is reached in 
about 20 min). Typical images of the storage area are 



FIG. 5: Images of a part of the storage area obtained under 
crossed polarizers. Droplets of adipic acid solution are stored 
in the channels. The concentration in adipic acid was grad- 
ually changed between the upper and the bottom channels. 
After crystallisation of all the droplets, a temperature gradi- 
ent is applied (low temperature on the left and high temper- 
ature on the right), (a) The dotted line separating droplets 
containing crystals from empty droplets give an estimation 
of the solubility limit, (b) Same image but with a different 
contrast displaying the droplets positions. 



presented on Fig. [H Since adipic acid crystals have bire- 
fringent properties, they are easily detected under crossed 
polarizers. The smallest detectable crystals size is about 
50x50 fim 2 at the magnification used. Figure E^a) en- 
ables us to directly observe the limit of crystal presence. 
Using interpolated temperature profiles such as the one 
displayed in Fig. [4j allows us to estimate the solubility 
temperatures for all the ten concentrations (we choose 
them in the middle of the two sucessive droplets with 
and without crystals). Figure [6] presents such solubility 
temperatures measured with our microfluidic device. The 
error corresponds to the temperatures difference between 
the two droplets enclosing the solubility limit positions. 
These results are in good agreements with data obtained 
from literature (27| . 

Naturally, the errors done on such measurements de- 
pend on the distance between two successive droplets, 
and on the amplitude of the temperature gradient. In 
our case, the temperature gradient of 0.7°C mm -1 and 
a typical distance of 3 mm between two droplets give an 
error of ±1°C. The application of smaller temperature 
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FIG. 6: (•) Solubility of adipic acid in water measured in the 
case of a temperature gradient of 0.7°C mm -1 , (o) Solubility 
data from literature, the dotted line is a guideline for eyes. 



gradients and the reduction of the distance between two 
successive droplets would give a better accuracy on the 
solubility limit. 

For the moment, the maximal temperature which can 
be investigated is limited by the evaporation of water 
through the PDMS layer [28|. Simple measurements 
show that the volume of an aqueous droplet stored in 
our device at 60° C, decreases by w 10% in 4 hr. Such an 
effect is negligible for the experiments described above 
(droplet filling time 20 min at 65° C), but may explain 
the small discrepancy observed on Fig. [6] at high temper- 



ature. We believe that the use of non-permeable materi- 
als such as glass, instead of PDMS, could easily broaden 
the possibilities offered by our system. 



IV. CONCLUSION 

In this work we have presented a new microfluidic 
tool to perform rapid screening of solubility diagrams. 
The device enables us to store hundreds of droplets 
(w 100 nL) of various chemical compositions in parallel 
microchannels, and to apply large temperature gradients. 
We have demonstrated using a model system (adipic acid 
in water), that we could easily and directly access to 
ten simultaneous measurements of the solubility curve 
on a large temperature range, in less than 1 hr, and with 
only 250 fiL of solution. To conclude, we believe our 
device is a suitable tool for solubility diagrams screen- 
ing, more rapid, with a better temperature control, and 
cheaper than classical robotic workstations. Such a mi- 
crofluidic tool may also be useful for many other applica- 
tions, where two-dimensional screening, temperature vs. 
composition, is required. 
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